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In-Depth Modeling 


• Mass balance (gas phase) 

- Solid to gas conversion (pyrolysis) 

- Gas advection (permeation) 

- Storage (in the internal pores) 

• Momentum balance 

- Balances friction with pressure 
gradient (Darcy’s model) 

• Energy balance 

- Conduction (in solid phase) 

- Storage (in solid and gas phases) 

- Pyrolysis energy 

- Advected energy (heat exchange 
with pyrolysis gases) 
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In-Depth Modeling 


• Pyrolysis 
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Arrhenius model for component 
extent-of-reaction 


Mass/Momentum Equation 
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- Neglected storage (quasi-steady) 
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Solution of M/M 
equation primarily 
needed for solution 
of energy equation. 
Analyst can choose 
level of fidelity. 


- 1-D simplification (with neglected storage) 
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Thermochemistry 
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ITRAC 
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• 1-D (planar, cylindrical, and spherical) 

• Variable-grid finite-volume method 

• Heat transfer, material pyrolysis, pore pressure, thermochemical ablation 

• Various mechanical erosion models 

• User-defined dynamic link libraries (DLLs) 

• Ignition model 



ITRAC vl.O 

Insulation Thermal Response and Ablation Code 
ATK Space Systems © 2008 
Promontory, Utah 
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Hero (Heat Transfer and Erosion Analysis Program) 
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• Multi-D 

• Variable-grid finite-element (hierarchical) method 

• Heat transfer, material pyrolysis, pore pressure, 
thermochemical ablation, raaiosity, structural 

• Fluid-thermal-structural interaction (FTSI) 
capabilities (FEM Builder) 

• Adaptive refinement 

• Parallel processing 





Hero Adaptive Refinement 



FTI Nozzle Analysis 


Run-Time Improvement with Parallel Processing 
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Hero (Heat Transfer and Erosion Analysis Program) 
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Modeling of Complex Ablation Scenario Nozzle Thermal Analysis Results 



Radiosity with Solar Effects 
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Material Properties for In-depth Models 

A premier aerospace and defense company 


Pyrolysis Properties 
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Arrhenius model for component 
extent-of-reaction 


Insulator Properties 

•Pyrolysis kinetics parameters 

•Virgin and char densities 

Pv Pc 
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•Virgin and char specific heat 
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Input for Surface Modeling 
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Boundary Conditions 
•Transport coefficient 
•Enthalpy (recovery and frozen) 
•Incident radiation heat flux 


Propellant Properties 
•Stanton number ratio 
•Elemental composition 


Surface Product Properties 

•Elemental composition 
•Enthalpy of products 


Insulator Properties 
•Char enthalpy (specific heat) 
and pyrolysis gas enthalpy 

•Radiation properties 
(emissivity and absorptivity) 
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Property Summary 


Pyrolysis kinetic parameters 
Density 

• Virgin and char 
Specific heat 

• Virgin and char 

• Versus temperature 
Thermal conductivity 

• Virgin and char 

• Versus temperature 
Elemental compositions 

• Virgin and char 

• Pyrolysis gas calculated 
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Heats-of-formation 

• Virgin and char 

• Pyrolysis gas value calculated 

• Used to calculate heat-of-pyrolysis 

• Used in surface thermochemistry 
Radiation properties 

• Emissivity and absorptivity 
For pore pressure 

• Porosity 

• Permeability 

• Pyrolysis gas molecular weight 

• Pyrolysis gas viscosity 



Structural Modeling 
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Structural modeling is often required for accurate assessment of structural integrity 
(including phenomena such as pocketing, ply-lifting, wedge-outs, and delamination) 

• Phenomena can have significant impact on thermal protection 

• Integrally tied to thermal responses 

• Accurate modeling of heat transfer, material pyrolysis, pore pressure, thermochemical 
ablation, etc. is critical 
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Structural - High Temperature Issues 
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Pocketing involves the expulsion of carbon-cloth phenolic (CCP) material due to a 
combination of stresses from thermal expansion and pore pressure driven stresses 

• Can significantly affect thermal protection 




Structural - High Temperature Issues 
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Ply lifting involves the lifting of a char cap of CCP material. Initial failure is caused 
by thermal expansion/contraction, lifting is caused by pore pressures 
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Structural - High Temperature Issues 
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Ply lifting char caps can slough, causing loss of thermal protection 

^371 03 sec I SRTM NC-6 24" Motor Nozzle Real-Time Radiography, 9/14/2006 1 

Pre-fire Marshall Space Flight Center NDE Team/EM20 



Structural Dependent on Thermal 
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Structural modeling is dependent on thermal modeling, for example: 

• Stiffness can be dependent on degree-of-char and presence of moisture 
- Charring causes material changes that affect the magnitude of the modulus 


- Moisture causes a plasticization of the matrix material 
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Strain 



Structural Dependent on Thermal 



Structural modeling is dependent on thermal modeling, for example: 

• Structural behavior is also influenced by pressures 
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Property Summary 


The following is a simple first principles equation that may be used for the structural 
model (complexity is usually added to address nonlinearities) 


= S m aZ + tJ„AP + a,, AT + ft, AM + ZlJ Av c 

1 .. H u M u M .. l\ , I 

y y y y y 


stress pore thermal moisture decomposition 

loading pressures expansion swelling 

Material properties summary (dependent on thermal state) 

Stiffness matrix (includes moduli and Poisson’s ratios) S,j k/ 

Pressure stress coupling, 77^ 

Coefficient of thermal expansion, a Y] 

Moisture expansion coefficient, 

Decomposition coefficient, x\\ 

Failure criteria 


Conjugate models are required for accurate simulation of thermal and structural behavior 
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On-going Work 
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Thermal 

• 3-D surface ablation in Hero 

• General thermal contact in Hero 

• Advanced surface thermochemistry code (ACE replacement) 

- Alumina impingement 

- Chemical kinetics 

• Material properties 

• Comprehensive validations against historical nozzle data 

Flow modeling 

Advanced carbon-cloth phenolic CCP modeling 

• Moisture 

• Structural coupling 
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